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Abstract 

We consider generalized linear stochastic dynamical systems with 
second-order state transition matrices. The entries of the matrix are 
assumed to be either independent and exponentially distributed or 
equal to zero. We give an overview of new results on evaluation of 
asymptotic growth rate of the system state vector, which is called the 
Lyapunov exponent of the system. 

1 Introduction 

The evolution of actual systems that occur in management, engineering, 
computer sciences, and other areas can frequently be represented through 
stochastic dynamic equations of the form 

z{k) = A{k)z{k - 1), 

where A{k) is a random state transition matrix, z{k) is a system state 
vector, and matrix- vector multiplication is thought of as defined in terms of 
a semiring with the operations of taking maximum and addition [HOIS]. 

In many cases, the analysis of a system involves evaluation of asymptotic 
growth rate of the system state vector z(k), which is normally referred to 
as the Lyapunov exponent [4;, 5]. 

Evaluation of the Lyapunov exponent typically appears to be a difficult 
problem even for quite simple systems. Related results include the solu- 
tions obtained in [6l [5] for systems with matrices of the second order with 
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independent and exponentially distributed entries. In [6], the Lyapunov ex- 
ponent is obtained in the case that all entries of the matrix are identically 
distributed with unit mean. 

Further results are given in [5j under the condition that the diagonal 
entries have one common distribution, whereas the ofF-diagonal entries do 
follow another distribution. A system with a matrix such that its diagonal 
entries are distributed with unit mean, and the off-diagonal entries are equal 
to zero is also examined. 

The purpose of this paper is to give an overview of new results which are 
related to evaluation of the Lyapunov exponent in generalized linear systems 
that have matrices of the second order with exponentially distributed entries 
(second-order exponential systems). 

2 Stochastic Linear Dynamical System 

Consider a dynamical system that can be represented through the linear 
equation in the semiring with the operations of maximum and addition 

z{k) = A{k)z{k - 1), 

where 




With ordinary notation, the above vector equation can be written as 

x{k) = max{x{k — 1) + a^, y{k — 1) + 
y{k) = max(x(A; - 1) +-fk,y{k - 1) + 6k). 

The Lyapunov exponent for the system is given by 

A= lim — max{x{k),y{k)). 

fc— i>oo k 

Suppose the sequences {a^}, {/3k}, {Tfc}, and {6k} each involve inde- 
pendent and identically distributed random variables; ak, A, 7m , and 6n 
are independent for any k, l,m,n. Finally, we assume that , (3k, jk, and 
6k have the exponential probability distributions with respective parameters 
H, V, cr, and r. 

3 Systems With Matrices Having Zero Entries 

We start with results obtained in [71 [8j for systems with state transition 
matrices having one or two nonrandom entries that are equal to zero. The 
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reduced number of random entries in the matrices allows one to simplify 
the evaluation of the Lyapunov exponent and usually gives quite compact 
results. 

The solution method is based on construction of a sequence of probability 
distribution functions. The convergence of the sequence is examined and the 
limiting distribution is derived as the solution of an integral equation. The 
Lyapunov exponent is then evaluated as the expected value of a random 
variable determined through the limiting distribution function. 

3.1 Matrix With Zero Off-Diagonal Entries 

The system state transition matrix together with its related result take the 
form 

With T = /i = 1 we have the result A = 1.25 which coincides with that 
in [5]. 

3.2 Matrix With Zero Diagonal 

In this case, the matrix and the Lyapunov exponent are represented as 
A{k) = ( ^,0, ,_4.^ + 7.. + 4.2 



Ik y ' Qua^u + a) 

3.3 Matrix With Zero Row or Column 

Provided that the second row in the matrix has only zero entries, we arrive 
at 



When the entries of the second column are zero, we have 



7fc 7 ' /icr(/i + cr)2(3^ + 4o-) 

3.4 Matrix With Zero Entry on Diagonal 

Consider a system with the state transition matrix 



A{k) 



Oik h 
Ik 



Whereas evaluation of the Lyapunov exponent for this system in the 
general case leads to rather cumbersome algebraic manipulations, there are 
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two main particular cases which offer their related results in a relatively 
compact form. Under the condition that cr = /x, we have 

~ 2Aiz/(36/x4 + U7n^iy + 215n'^u^ + 130/xz/3 + 28z/4) " 
Provided that a = v, the value of the Lyapunov exponent is given by 

A = P{n,iy)/Q{i2,u), 

where 

P(/x, z/) = 15/ + + 624/x^i^2 + 1382^^1^^ + 1838/i/^ + 1592fi^u^ 

+ 973/x^z/^ + 384/xi/^ + 64I/^ 
Q(/x, z/) = nu{fi + z/)2(12/x^ + + 286/x^i/2 + 397/x2z/^ + 256/xz/'' + 64z/^). 

3.5 Matrix With Zero Entry Below Diagonal 

Suppose that there is a system with the state transition matrix defined as 



Aik) 



Oik h 
5k 



Consider two particular cases. With the condition u = fj,, we have the 
result 

A = P(/x,r)/Q(/x,r), 

where 

P(/x, r) = 288/ + 1048/r + 1936/t2 + 2688/xV^ + 3012/r'' 
+ 2226/r^ + 941/xV^ + 204/xr^ + 17r^, 

Q(/x, r) = 2/xr(144/ + 524/r + 968/r^ + 1200/t^ + 910/xV 
+ 387/xV^ + 84/xr^ + 7t^). 

Provided that t = ji, the solution takes the form 

A = P{ii,u)/Q{n,u), 

where 

Pill, u) = 256//^° + 2112/xV + 8044/i/2 + 19355/z/^ + 32167/z/^ 
+ 36887/z/^ + 28709/z/^ + 14854/z/^ + A912f?i^^ 
+ 4/xz/^ + 80z/^°, 

Q(/x, r^) = 2/xz/(/x + i^)(192/ + 1344/ z/ + 4047/ z^^ + 6770/z/^ 

+ 6799/z/'' + 4216/z^^ + 1600/z/^ + 344/xz/^ + 32z/^). 
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4 General Second-Order Exponential System 



Consider a general second-order exponential system which has the matrix 

^<^> ^ ( : t 

with its entries , (3k, 7fc , and 6k assumed to be independent random 
variables that are exponentially distributed with the respective parameters 
/i, v , a, and r. 

To get the Lyapunov exponent, a computational technique developed 
in [9j can be implemented which reduces the problem to the solution of 
a system of linear equations, accompanied by the evaluation of a linear 
functional of the system solution. The technique leans upon construction 
and examination of a sequence of probability density functions. It is shown 
that there is a one-to-one correspondence between the density functions and 
vectors in a vector space. The correspondence is then exploited to provide 
for the solution in terms of algebraic computations. 

Based on the above technique, the Lyapunov exponent can be evaluated 
as follows. First we introduce the vectors 



CJl — (a;io,Wll,Wl2,^13) ) ^^2 — (W20,W21,U;22, W23) 

Furthermore, we define the matrices 
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Suppose that the vector cj is the solution of the system 



(/ - W)u = 0, 

CJio + 1^20 = 1; 



where I represents identity matrix, and 

W = 



UiVu UiVi2 
U2V21 U2V22 



The value of the Lyapunov exponent is then given by 

A = qjui + qju;2, 



where qi and 92 are vectors such that 



\ 



T{fl+T)((T+T){fl+a+T) 

\ (u+T)(p+u+T){u+a+T)ifi+u+a+T) I 



( 



VT(y-\-T) 
f T(2/l+i/ + T) 



_ /i(/i+i^)(/j+r)(^i+!/+r) 

y2 — VT(y^2a^T) 

(T(!/+(T)(o-+r)(z^+(T+r) 

1 t^r(2M+i/+2o-+r) . 

\ (At+o-)(/i+i/+cr)(At+o-+T)(^j+i/+cr+r) / 

Suppose that t = [i and a = v . Implementation of the above technique 
gives the solution 



where 



160/i^° + 1776^^ + 8220/i/2 + 2\2n%pJ + SSSOS/i'^z^^ 



+ GDO/i^l/^ + 321^2j,6 _^ gQ^j^7 _^ gj^8)_ 

Note that the obtained solution coincides with that in [S]. 
The author is grateful to the anonymous reviewer for valuable comments 
and suggestions. 
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